Non-symmetric trapped surfaces in the Schwarzschild and Vaidya spacetimes

Marginally trapped surfaces (MTSs) are commonly used in numerical relativity to locate black holes. For dynamical black holes, it is not known generally if this procedure is sufficiently reliable. Even for Schwarzschild black holes, Wald and Iyer constructed foliations which come arbitrarily close to the singularity but do not contain any MTSs. In this paper, we review the Wald-Iyer construction, discuss some implications for numerical relativity, and generalize to the well known Vaidya spacetime describing spherically symmetric collapse of null dust. In the Vaidya spacetime, we numerically locate non-spherically symmetric trapped surfaces which extend outside the standard spherically symmetric trapping horizon. This shows that MTSs are common in this spacetime and that the event horizon is the most likely candidate for the boundary of the trapped region.Comment: 4 pages, 3 figures; v2: minor modifications; v3: clarified conclusion

Automatic correction of hand pointing in stereoscopic depth

In order to examine whether stereoscopic depth information could drive fast automatic correction of hand pointing, an experiment was designed in a 3D visual environment in which participants were asked to point to a target at different stereoscopic depths as quickly and accurately as possible within a limited time window (≤300 ms). The experiment consisted of two tasks: "depthGO" in which participants were asked to point to the new target position if the target jumped, and "depthSTOP" in which participants were instructed to abort their ongoing movements after the target jumped. The depth jump was designed to occur in 20% of the trials in both tasks. Results showed that fast automatic correction of hand movements could be driven by stereoscopic depth to occur in as early as 190 ms.This work was supported by the Grants from the National Natural Science Foundation of China (60970062 and 61173116) and the Doctoral Fund of Ministry of Education of China (20110072110014)

High accuracy binary black hole simulations with an extended wave zone

We present results from a new code for binary black hole evolutions using the moving-puncture approach, implementing finite differences in generalised coordinates, and allowing the spacetime to be covered with multiple communicating non-singular coordinate patches. Here we consider a regular Cartesian near zone, with adapted spherical grids covering the wave zone. The efficiencies resulting from the use of adapted coordinates allow us to maintain sufficient grid resolution to an artificial outer boundary location which is causally disconnected from the measurement. For the well-studied test-case of the inspiral of an equal-mass non-spinning binary (evolved for more than 8 orbits before merger), we determine the phase and amplitude to numerical accuracies better than 0.010% and 0.090% during inspiral, respectively, and 0.003% and 0.153% during merger. The waveforms, including the resolved higher harmonics, are convergent and can be consistently extrapolated to $r\to\infty$ throughout the simulation, including the merger and ringdown. Ringdown frequencies for these modes (to $(\ell,m)=(6,6)$) match perturbative calculations to within 0.01%, providing a strong confirmation that the remnant settles to a Kerr black hole with irreducible mass $M_{\rm irr} = 0.884355\pm20\times10^{-6}$ and spin $S_f/M_f^2 = 0.686923 \pm 10\times10^{-6}

Multi-patch methods in general relativistic astrophysics - I. Hydrodynamical flows on fixed backgrounds

Many systems of interest in general relativistic astrophysics, including neutron stars, accreting compact objects in X-ray binaries and active galactic nuclei, core collapse, and collapsars, are assumed to be approximately spherically symmetric or axisymmetric. In Newtonian or fixed-background relativistic approximations it is common practice to use spherical polar coordinates for computational grids; however, these coordinates have singularities and are difficult to use in fully relativistic models. We present, in this series of papers, a numerical technique which is able to use effectively spherical grids by employing multiple patches. We provide detailed instructions on how to implement such a scheme, and present a number of code tests for the fixed background case, including an accretion torus around a black hole.Comment: 26 pages, 20 figures. A high-resolution version is available at http://www.cct.lsu.edu/~bzink/papers/multipatch_1.pd

The Effect of Pictorial Illusion on Prehension and Perception’,

Abstract The present study examined the effect of a size-contrast illusion (Ebbinghaus or Titchener Circles Illusion) on visual perception and the visual control of grasping movements. Seventeen right-handed participants picked up and, on other trials, estimated the size of "poker-chip" disks, which functioned as the target circles in a three-dimensional version of the illusion. In the estimation condition, subjects indicated how big they thought the target was by separating their thumb and fore nger to match the target's size. After initial viewing, no visual feedback from the hand or the target was available. Scaling of grip aperture was found to be strongly correlated with the physical size of the disks, while manual estimations of disk size were biased in the direction of the illusion. Evidently, grip aperture is calibrated to the true size of an object, even when perception of object size is distorted by a pictorial illusion, a result that is consistent with recent suggestions that visually guided prehension and visual perception are mediated by separate visual pathways

Hyperboloidal slices for the wave equation of Kerr-Schild metrics and numerical applications

We present new results from two open source codes, using finite differencing and pseudo-spectral methods for the wave equations in (3+1) dimensions. We use a hyperboloidal transformation which allows direct access to null infinity and simplifies the control over characteristic speeds on Kerr-Schild backgrounds. We show that this method is ideal for attaching hyperboloidal slices or for adapting the numerical resolution in certain spacetime regions. As an example application, we study late-time Kerr tails of sub-dominant modes and obtain new insight into the splitting of decay rates. The involved conformal wave equation is freed of formally singular terms whose numerical evaluation might be problematically close to future null infinity.Comment: 15 pages, 12 figure

Accurate evolutions of inspiralling neutron-star binaries: prompt and delayed collapse to black hole

Binary neutron-star (BNS) systems represent primary sources for the gravitational-wave (GW) detectors. We present a systematic investigation in full GR of the dynamics and GW emission from BNS which inspiral and merge, producing a black hole (BH) surrounded by a torus. Our results represent the state of the art from several points of view: (i) We use HRSC methods for the hydrodynamics equations and high-order finite-differencing techniques for the Einstein equations; (ii) We employ AMR techniques with "moving boxes"; (iii) We use as initial data BNSs in irrotational quasi-circular orbits; (iv) We exploit the isolated-horizon formalism to measure the properties of the BHs produced in the merger; (v) Finally, we use two approaches, based either on gauge-invariant perturbations or on Weyl scalars, to calculate the GWs. These techniques allow us to perform accurate evolutions on timescales never reported before (ie ~30 ms) and to provide the first complete description of the inspiral and merger of a BNS leading to the prompt or delayed formation of a BH and to its ringdown. We consider either a polytropic or an ideal fluid EOS and show that already with this idealized EOSs a very interesting phenomenology emerges. In particular, we show that while high-mass binaries lead to the prompt formation of a rapidly rotating BH surrounded by a dense torus, lower-mass binaries give rise to a differentially rotating NS, which undergoes large oscillations and emits large amounts of GWs. Eventually, also the NS collapses to a rotating BH surrounded by a torus. Finally, we also show that the use of a non-isentropic EOS leads to significantly different evolutions, giving rise to a delayed collapse also with high-mass binaries, as well as to a more intense emission of GWs and to a geometrically thicker torus.Comment: 35 pages, 29 figures, corrected few typos to match the published version. High-resolution figures and animations can be found at http://numrel.aei.mpg.de/Visualisations/Archive/BinaryNeutronStars/Relativistic_Meudon/index.htm

Introduction to dynamical horizons in numerical relativity

This paper presents a quasi-local method of studying the physics of dynamical black holes in numerical simulations. This is done within the dynamical horizon framework, which extends the earlier work on isolated horizons to time-dependent situations. In particular: (i) We locate various kinds of marginal surfaces and study their time evolution. An important ingredient is the calculation of the signature of the horizon, which can be either spacelike, timelike, or null. (ii) We generalize the calculation of the black hole mass and angular momentum, which were previously defined for axisymmetric isolated horizons to dynamical situations. (iii) We calculate the source multipole moments of the black hole which can be used to verify that the black hole settles down to a Kerr solution. (iv) We also study the fluxes of energy crossing the horizon, which describes how a black hole grows as it accretes matter and/or radiation. We describe our numerical implementation of these concepts and apply them to three specific test cases, namely, the axisymmetric head-on collision of two black holes, the axisymmetric collapse of a neutron star, and a non-axisymmetric black hole collision with non-zero initial orbital angular momentum.Comment: 20 pages, 16 figures, revtex4. Several smaller changes, some didactic content shortene